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tJ.  ABSTRACT 

The  differential  and  total  cross  sections  for  ellipsoids  and  elliptic  cylinders 
having  Gaussian  electron  density  distributions  have  been  obtained  ty  a  ray  trac¬ 
ing  procedure.  Calculations  for  the  case  of  an  external  magnetic  field  were  re-  f 
stricted  to  the  ellipsoidal  distributions.  The  results  show  that  the  scattering  is 
extremely  sensitive  to  the  orientation  of  the  body.  A  peak  in  cross  section  oc¬ 
curs  at  the  scattering  angle  corresponding  to  the  ray  normal  to  the  critical  sur¬ 
face,  and  increases  as  the  surface  becomes  flatter.  The  cross  section  is  sensi¬ 
tive  to  the  ratio  of  peak  density  to  critical  density  for  moderate  values  but  be¬ 
comes  relatively  insensitive  when  the  ratio  exceeds  3.  The  total  cross  section 
is  a  very  sensitive  function  of  both  orientation  and  ratio  of  major  to  minor  axes. 
The  introduction  of  a  magnetic  field  decreases  the  ordinary  ra1  cross  section; 
the  extraordinary  ray  exhibits  higher  values  only  in  the  forwai  d  sc'.ctering  re¬ 
gion,  but  is  always  higher  for  the  spherical  case.  Comparison  of  the  Gaussian 
ellipsoid  with  the  corresponding  conducting  ellipsoid  shows  that  the  Gaussian 
has  a  large  cross  section  in  the  forward  region  but  considerably  lower  values  in 
the  backscatter  region. 
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Abstract 

The  differential  and  total  cross  sections  for  ellipsoids  and  elliptic  cylinders 
having  Gaussian  electron  density  distributions  have  been  obtained  by  a  ray  tracing 
procedure.  Calculations  for  the  case  of  an  external  magnetic  field  were  restricted 
to  the  ellipsoidal  distributions.  The  results  show  that  the  scattering  is  extremely 
sensitive  to  the  orientation  of  the  body.  A  peak  in  cross  section  occurs  at  the  scat¬ 
tering  angle  corresponding  to  the  ray  normal  to  the  critical  surface,  and  increases 
as  the  surface  becomes  flatter.  The  cross  section  is  sensitive  to  the  ratio  of  peak 
density  to  critical  density  for  moderate  values  but  becomes  relatively  insensitive 
when  the  ratio  exceeds  3.  The  total  cross  section  is  a  very  sensitive  function  of 
both  orientation  and  ratio  of  major  to  minor  axes.  The  introduction  of  a  magnetic 
field  decreases  the  ordinary  ray  cross  section;  the  extraordinary  ray  exhibits 
higher  values  only  in  the  forward  scattering  region,  but  is  a’ways  higher  for  the 
spherical  case.  Comparison  of  the  Gaussian  ellipsoid  with  the  corresponding  con¬ 
ducting  ellipsoid  shows  that  the  Gaussian  has  a  large  cross  section  in  the  forward 
region  but  considerably  lower  values  ;n  the  backscatter  region. 
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Scattering  of  HF  Radio  Waves  by  Elliptical 
Electron  DenHty  Distributions 


1.  INTRODUCTION 

In  recent  years  there  has  been  ocniG  work  performed  on  the  scattering  of  High 
Frequency  (HF)  radio  waves  by  artificial  charge  distributions  in  the  ionosphere. 
The  most  notable  of  these  are  the  "SECEDE"  Bat  ium  Releases  which  have  been 
discussed  in  the  literature  [Bates,  1971;  Rao,  et  al.,  1971;  T  iome,  196&),  Up  to  this 
time,  there  has  not  been  a  thorough  investigation  into  the  scattering  hy  these  re¬ 
leases  as  a  function  of  angle  of  the  incident  wave  to  the  major  axis  of  the  releases. 

We  have  therefore  calculated  the  total  cross  section  t  nd  the  differential  cross 
section  ot  a  long  cylindrical  charge  d.stributim  with  an  ellipsoidal  cro3a  section 
and  a  two-dimensional  Gaussian  electron  denary  distribution  for  several  orienta¬ 
tions  of  the  m-jor  axis  of  the  distribution  to  the  incident  wave.  To  uerform  these 
calculate  a  w*.  apply  the  theory  of  ray  optics  am  use  I.aa^lgro  'e's  differential 
equations  (Kelso,  1964)  to  calculate  the  ray  paths.  To  calculate  the  index  ot  refrac¬ 
tion,  we  have  employed  the  Appleton- Hartree  dispersion  equation  i  Kelso,  1964) 
neglecting  absorption  by  the  medium. 

In  the  following  discussion,  v/e  will  present  cur  results  for  -he  differential  and 
total  cross  section  of  long  cylindrical  bodies  wnich  have  various  ellipsoic.al  cross 
sections  and  differing  electron  contents.  In  addition,  we  will  treat  the  case  of  a 
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prolate  ellipsoid  when  the  direction  of  propagation  of  the  incident  wave  is  along  the 
major  axis  of  the  ellipsoid. 

Although  it  is  known  that  the  artificially  induced  charge  distributions  take  an 
ellipsoidal  form,  we  will  initially  treat  the  simpler  case  of  an  infinitely  long  cylin¬ 
der  with  an  elliptic  cross  section.  The  special  case  of  the  incident  plane  wave  par¬ 
allel  to  the  major  axis  of  the  charge  distribution  can  be  treated  as  a  prolate  ellip¬ 
soid  since  there  is  no  dependence  on  azimuthal  angle,  and  thus  the  cros^  section  is 
a  function  of  scattering  angle  only.  The  azimuthal  dependence  of  the  cross  section, 
which  is  important  in  the  case  of  an  incident  wave  at  an  oblique  angle  to  the  charge 
distribution,  has  not  been  included  in  our  calculations. 


2.  ANALYSIS 

For  a  spherically  symmetric  distribution  the  scattering  cross  section  is  given 
by  the  formula  (Merzbacher,  1961) 

■  SETT  l-SH  «> 

where  b  is  the  impact  parameter,  and  r  he  scattering  angle.  We  have  introduced 
a  factor  of  4ff  in  order  to  be  consistent  with  the  definition  used  in  radar  cross  sec¬ 
tion  studies. 

A  more  general  formula  (Merzbacher,  1961)  for  the  differential  cross  section  is 

dA  =  a(«)dA/  (2) 

where  dA  is  the  incident  flux  per  unit  area  and  dA7  the  area  into  which  it  is 
scattered. 

For  a  long  cylinder,  neglecting  end  effects, 

dA  =  db  dy  (3) 

where  b  is  the  impact  parameter  and  y  is  a  length  along  the  cylinder.  In  cylin¬ 
drical  coordinates,  we  have 

o(U)  dA*  =  c(0,y)  d0  dy  (4) 

where  6  ;s  the  scattering  angle.  Substituting  Eqs.  (3)  and  (4)  into  Eq.  (2), 


dbdy  =  0(6, y)  d0dy. 


(5) 
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Since  the  cylinder  has  a  uniform  cross  sectional  area  along  its  length,  0  is  a 
function  only  of  6,  and  we  can  thus  integrate  the  y  component  to  obtain 

db  =  0(9  )d0.  (8) 

Introducing  the  factor  of  2ir  to  be  consistent  with  radar  cross  sections,  we 
obtain 

0(B)  =  2ir  |-^-|.  (7) 

Eq.  (7)  was  used  to  calculate  the  differential  cross  section  per  unit  length  of 
the  cylinder.  For  convenience,  we  have  taken  the  unit  of  length  to  be  1  km.  The 
total  cross  section  can  be  calculated  by  integrating  the  above  equation  over  6. 

In  order  to  evaluate  the  derivative  in  Eq.  (1 ),  we  must  calculate  the  paths  of  a 
number  of  rays  in  order  to  obtain  the  dependence  of  b  on  6.  These  ray  paths  are 
determined  by  Haselgrove's  equations,  which  we  will  discuss  in  the  next  section. 
The  evaluation  of  the  derivative  was  worked  out  with  the  help  of  Rosenberg  (Rosen¬ 
berg,  1971),  who  has  written  a  spline  fitting  program  (Ahlberg,  1967)  which  allows 
us  to  fit  a  function  to  a  number  of  cubic  polynomials  so  that  the  first  derivative  is 
continuous  or  nearly  continuous  at  each  point  along  the  curve.  Thus  we  can  calcu¬ 
late  o(6)  under  a  variety  of  different  conditions. 


3.  PROCEDURE 

We  now  discuss  the  evaluation  of  Haselgrove's  equations  and  the  computer  code 
which  was  written  to  calculate  the  scatteiing  cross  section  of  a'long  cylinder  with 
an  ellipsoidal  cross  section.  Figure  1  shows  the  system  of  coordinates  which  we 
have  used  in  our  analysis.  The  x-axis  is  horizontal,  the  z-axis  is  in  the  vertical 
plane,  and  the  y-axis  points  into  the  paper.  The  rotation  angle  of  the  ellipse  0^  is 
measured  counterclockwisr  from  the  positive  x-axis  to  the  major  axis  of  the  el¬ 
lipse.  The  impact  parameter  b  is  measured  from  the  x-axis  and  is  incident  from 
the  left.  The  perimeter  of  the  ellipse  is  the  curve  along  which,  when  no  magnetic 
field  is  present,  the  index  of  refraction  ji(x,z)  is  zero  (critical  ellipse).  The  dis¬ 
tance  Z  is  the  z-coordinate  of  the  point  in  the  left-hand  plane,  where  the  tangent  to 
the  ellipse  is  perpendicular  to  the  x-axis. 

In  order  to  solve  F,q.  (7),  a  Fortran  computer  program  was  written  to  calculate 
b(0)  and  subsequently  evaluate  <3(9)  for  a  number  of  different  parameters.  Follow¬ 
ing  is  a  brief  description  of  the  steps  involved  in  the  computations.  Appendix  A 
contains  a  copy  of  the  computer  code  which  was  written. 
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The  first  step  in  the  calculation 
of  cross  section  is  the  determination 
of  the  rays  which  define  b(0).  In 
order  to  calculats  the  rays,  we  need 
a  set  of  impact  parameters  which 
should  cover  the  complete  range  of 
scattering  angles.  This  set  of  impact 
parameters  varies  according  to  the 
orientation  and  width  of  the  cylinder. 
To  determine  a  complete  set  of  im¬ 
pact  parameters,  we  calculate  a  quan¬ 
tity  t!  about  which  we  symmetrically 
distribute  the  impact  parameters  by 

Z  =  Z  -  .3  sin(0pj)  2<niax  (B) 

where  is  the  angle  the  major  axis 
makes  with  the  x-axis,  and  Zmax  is  the  maximum  z  coordinate  of  the  critical  el¬ 
lipse.  The  above  empirical  formula  was  found  to  work  sufficiently  well  for  the  ro¬ 
tation  angles  0°,  45°  and  90°. 

The  impact  parameters  were  closely  spaced  near  the  value  of  Z#,  but  the  spac¬ 
ing  increases  as  we  proceeded  away  from  the  line  of  symmetry.  Tne  maximum 
sep?  ation  between  impact  r  >.rameters  was  .25  km,  and  they  extended  to  a  distance 
of  4  Zmax  from  Z. 

Once  we  have  chosen  the  impact  parameters,  we  launch  each  ray  from  an  x 
coordinate  of  -5  km  towards  tne  ellipse.  Haselgrove’s  equations,  which  define  the 
path  that  a  ray  will  take  m  a  medium  with  an  index  of  refraction  p  are: 

“ ■  =  c<>s  a  +  sin  oi) 

dt  if  da 

—  =  -%  (p  sin  a  -  sin  a)  (9) 

dt  p4  da 

-  ~  (cos  a  -  sin  a  — ) 
dt  p^  dz  dx 

where  p(x,z)  is  the  index  of  refraction  at  the  point  (x,z/.  The  phase  angle  of  the  ray 

Az 

is  a ,  and  is  the  angle  between  in?  ray  direction  and  the  initial  direction  of 


Figure  1.  Geometry  and  Coordinate  Sys¬ 
tem  of  the  Electron  Density  Distribution 
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propagation.  To  solve  Haselgrove's  equations  by  computer  techniques,  Eqs.  (9) 
were  transformed  into  Eqs.  (10), 

Ax  =  Ad  0*  cos  a  +  ||  sin  a)  /  fi2 

Az  =  Ad  ( n  sin  a  -  cos  a)  /  a2  (10) 

Aa  =  Ad  (cos  a  -  sin  a  |£)  /  \f 


where  Ad  =  c  At  and  c  is  the  speed  of  light.  The  quantity  Ad  which  was  used  in 
integrating  the  above  equations  represents  the  step  size  which  at  a  maximum  was 
0.5  km.  This  value  was  decreased  by  the  program  such  that  in  no  step  would  the 
phase  angle  change  by  more  than  3°.  The  integration  was  carried  out  by  means  of 
a  Rungc-Kutta  method  for  the  solution  of  a  set  of  simultaneous  linear  differential 
equations  (Scarborough,  1930).  When  the  calculation  of  the  rays  was  performed 
with  a  step  size  of  0.25  km,  half  the  normal  step  3ize,  the  maximum  difference  in 
the  scattering  angle  was  0.009°.  The  difference  in  the  scattering  cross  section  was 
less  than  0.001  km2  where  the  maximum  scattering  error  occurred. 

The  value  of  p  at  each  point  that  the  equations  were  evaluated  is  given  by  the 
Appleton- Hartr>. 2  Dispersion  Equation, 


(x,z )  =  1  - 


X(l-X» 


1  -  X  -  ±  j Y2  +  (1  -  X)2 

TV  T 


“here 


X(x,z)  =  p  exp 


'  X2  z2  > 

T  *~T 

x„  z„  ; 

o  o  / 


(11) 


(12) 


A  measure  of  the  hardness  of  the  charge  distribution  is  p,  which  is  defined  as 
the  ratio  of  peak  plasma  density  to  the  critical  plasma  density.  If  p  is  somewhat 
greater  than  1,  we  have  a  hard  charge  distribution,  for  which  we  consequently  have 
a  significant  region  of  backscatter. 

The  Gaussian  half  widths  of  the  charge  distributions  are  x  and  z  . 

oo 

The  gyromagnetic  ratio  (ratio  of  ion  frequency/ incident  frequency)  is  Y,  and 
Yj  and  Y^  are  che  transverse  and  longitudinal  components  of  the  gyromagnetic 
ratio, 
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Yt  =  Y  sin  (a  -  $H) 
Yj  =  Y  cos  (a  -  flH) 


(13) 


where  0jj  is  the  angle  the  magnetic  field  makes  with  the  incident  wave.  In  the  case 
of  a  nonmagnetic  field,  Y  =  0,  and  the  equation  for  the  index  of  refraction  sim¬ 
plifies  to 

M2(x,z)  =  1  -  X.  (14) 

Once  each  impact  parameter  has  traced  a  ray,  the  resulting  table  of  b  and  9  define 
b  (9). 

The  value  of  b  for  9  =  180°  is  defined  as  the  displacement,  and  it  is  subtracted 
from  all  values  of  b.  The  reason  for  the  displacement  is  the  asymmetry  of  an  el¬ 
liptical  charge  distribution  when  it  is  rotated  with  respect  to  the  incident  plane 
wave,  or  when  a  magnetic  field  is  present. 

The  table  of  b  and  9  is  then  interpolated  at  intervals  of  1  °  with  a  cubic  spline, 
and  the  slope  at  each  point  is  calculated.  Since  the  slope  has  a  few  discontinuities 
due  to  imperfections  in  the  method  of  interpolation,  it  is  smoothed  by  taking  an  11 
point  running  average  twice. 

The  process  of  smoothing  is  of  questionable  value  since  the  total  cross  sec¬ 
tions  appear  to  be  increased  by  as  much  as  1 0%  when  the  smooth  data  is  used.  The 
effect  of  smoothing  appears  to  "lift"  the  curve,  thus  increasing  the  value  of  the  in¬ 
tegral.  The  region  where  the  smoothing  is  needed  is  normally  in  the  backscatter 
region  where  the  cross  section  is  small.  Thus  it  would  seem  that,  although  not 
aesthetically  pleasing,  the  curves  with  small  discontinuities  in  the  backscatter  re¬ 
gion  are  more  accurate  over  the  completed  range  of  angles  than  the  smoothed 
curves. 

Once  we  have  the  derivative,  the  cross  sect. on  is  easily  calculated  by  Eq.  (7), 
and  the  total  cross  section 


The  values  at  0°  and  360°  are  determined  by  linear  extrapolation. 

In  order  to  compare  the  cross  section  values  to  a  cylindrical  rod  of  circular 

cross  section,  we  have  calculated  the  values  of  x  and  z  such  that  the  total  elec- 
’  o  o 

tron  content  in  a  slice  of  the  cylinder  is  the  same. 
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The  electron  content  in  a  circular  cylinder  of  unit  length  is 


dx  dz  =  ptr 


(16) 


and  the  electron  content  in  an  ellipsoidal  cylinder  of  unit  length  is 


pv. 


■>zo ' 


(17) 


Therefore  the  electron  content  for  a  slice  of  a  cylinder  will  remain  constant  if  the 

product  of  the  major  and  minor  axes  remains  constant. 

If  w<-  take  r  =1  km,  then 
o  1 

1  =  x  z  .  (18) 

o  o 

If  we  want  an  ellipse  with  x  /z  =  5  =  R  ,  then 

0  0  o 

1  =  5z  2.  (19) 

o 

We  will  now  be  able  to  determine  not  only  how  the  orientation  of  the  charge  distri¬ 
bution  changes  the  cross  sections,  but  how  the  shape  of  the  distribution  affects  the 
cross  section. 

The  program  which  was  written  had  a  number  of  parameters  which  could  be 
changed  to  investigate  their  effects  on  the  cross  section.  These  parameters  were 
p,  xq,  zQ  and  6^.  The  output  of  the  program  was:;  impact  parameters,  the  total 
cross  section  per  unit  length,  and  a  tabulation  of  b,  db/d0,  o(6)  at  1°  intervals,  and 
total  cross  section. 


4.  RESULTS  FOR  A  LONG  CYLINDER 


The  program  to  calculate  the  cross  sections  produces  6  graphs  for  each  choice 
of  parameters.  The  curves  are  the  cay  paths  (Figures  2a,  3a,  4a),  the  function  f (0) 
(Figures  2b,  3b,  4b),  the  derivative  db/d0  (Figures  2c,  3c,  4c),  and  the  cross  sec¬ 
tion  (Figures  2d,  3d,  4d).  In  addition,  the  smooth  values  of  derivative  and  cross 
section  are  also  plotted  but  are  not  shown.  Figure  2  has  the  ellipse  unrotated. 
Figure  3  has  the  ellipse  rotated  45°,  and  Figure  4  has  the  ellipse  rotated  9*/  >-i 

all  three  cases,  p  =  2  and  Rq  =  3.  In  Figures  5,  6,  7  we  have  shown  the  differen¬ 
tial  cross  section  for  R0  =  1,  2,  and  5;  Or  =  0,  45°  and  90°  with  p  -  2.0. 


Figure  2.  Scattering  From  an  Elliptic  Cylinder  with  p  =  2.0,  R  =3,  and  =  0  ; 
(a)  Ray  Paths,  (b)  b  vs  6,  (c)  db/d0,  (d)  o(9)  0  K 


From  the  figures,  it  can  be  seen  how  the  detail  structure  varies  with  different 
orientations,  and  the  sharp  peak  when  Rq  ^  largest. 

If  we  carefully  examine  the  differential  cross  section  curves  in  Figures  6  and  7, 
we  notice  that  the  cross  section  reaches  a  relative  maximum  between  the  end  points. 
This  peak,  which  is  dependent  on  the  angle  between  the  major  axis  of  the  ellipse  and 
the  propagat.on  direction,  increases  the  total  cross  section.  The  radius  of  curva¬ 
ture  of  the  charge  distribution  is  also  an  important  factor.  The  peak  is  highest 
when  the  curvature  is  3mall,  and  the  wave  normal  is  perpendicular  to  the  tangent 
of  the  charge  distribution  at  the  point  of  contact  with  the  critical  surface. 
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Figure  4.  Scattering  From  an  Elliptic  Cylinder  with  p  =  2.0,  R  =  3,  and  0R  =  90°; 
(a)  Ray  Path3,  (b)  b  vs  9,  (c)  db/dP,  (d)  o(6)  0  K 


Figure  5.  Differential  Cross  Section  for  p  =  2.0,  9„  =  0°,  and  Several  Values  of  R 
(a)  R  -  1,  (b)  R  2,  (c)  R  -  5  K  ' 


Figures.  Differential  Cross  Section  for  p  =  2.0,  =  45°,  and  Several  Values  of  R( 

(a)  RQ  =  1,  (b)  R  =  2,  (r)  R  =  5 


Figure  7.  Differential  Cross  Section  for  p  =  2.0,  9r>  =  90°, 
(a)  R0  =  1,  (b)  R0  =  2,  (c)  R0  =  5 
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Id  Figure  8  we  have  plotted  the  location  of  the  peak  vs  the  rotation  angle  of  the 
ellipses  for  Rq  =  2,  3  and  5.  The  points  which  do  not  lie  on  the  line  =  2 0R 
probably  are  different  due  to  the  reflecting  body  being  an  ellipsoid. 


Figure  8.  Angle  of  Relative  Maximum  of  Differ¬ 
ential  Cross  Section  vs  0^ 


In  Figure  9  we  have  shown  how  the  cross  section  at  an  orientation  0°,  45° 
and  90°  changes  as  we  vary  p  and  keep  Rq  =  3.  One  of  the  important  features 
we  have  noted  is  the  insensitivity  of  the  total  cross  section  to  p.  Figure  9  shows 
the  variation  in  cross  section  is  not  very  significant  for  p  greater  than  2. 

In  Figure  10  we  have  plotted  the  total  cross  section  per  unit  length  vs  Rq  for 
four  different  angles  of  rotation.  This  graph  demonstrates  how  elongation  of  the 
distribution  can  increase  or  decrease  the  total  cross  sect  .on;'  for  an  angle  of  about 
20°,  the  cross  section  remains  relatively  constant  regardless  of  the  value  or  Rq. 

In  Figure  11  we  can  see  the  large  change  in  the  total  cross  section  as  a  func¬ 
tion  of  angle.  This  graph  demonstrates  the  importance  of  knowing  how  the  ellipse 
is  aligned  in  order  to  properly  interpret  radar  data.  In  addition,  we  have  also 
plotted  the  smooth  data  on  this  curve  for  R  =  5  to  show  how  the  total  cross  sec¬ 
tion  has  increased  with  the  smoothing. 
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Figure  11.  Total  Cross  Section  vs  0R  for  p  =  2.0, 
and  Several  Values  of  R 

o 

5.  RESULTS  FOR  THE  PROLATE  ELLIPSOID 

In  the  case  of  a  prolate  ellipsoid,  the  cross  section  of  the  electron  density  dis¬ 
tribution  perpendicular  to  the  major  axis  is  a  circle.  If  we  have  a  plane  wave  in¬ 
cident  along  the  major  axis  of  the  ellipse,  it  encounters  a  circularly  symmetric 
charge  distribution,  and  we  can  apply  Eq.  (1 )  to  calculate  the  differential  cross 
section. 

In  Figure  12  we  show  the  differential  cros3  section  for  R  =1,2  and  5.  From 
these  curves  it  appears  that  elongating  the  distribution  has  an  effect  similar  to  de¬ 
creasing  the  size  of  a  spherical  distribution  or  decreasing  the  value  of  p.  If  we 

2  v 

1  n 

integrate  Eq.  (1),  we  obtain  the  total  cross  section  a  =  ^  \  tf(0)d0  of  the  ellipsoid. 

o 

In  Figure  13  we  have  plotted  the  total  cross  section  vs  R0  keeping  the  total 
electron  content  constant.  For  this  case  Eq.  (19)  has  to  be  modified  such  that  zQ  is 
a  function  of  the  cube  root  of  the  ratio.  From  Figure  13  we  can  see  the  dramatic 
decrease  in  the  total  cross  section  as  the  incoming  wave  encounters  a  smaller  area 
perpendicular  to  the  direction  of  propagation  and  thus  has  a  smaller  backscatter 
region. 


IV 


Figure  12.  Differential  Cross  Section  vs  Scattering  Angle  with  Azimuthal  Symmetry, 

p  =  2.0,  0R  =  0°  and  Several  Values  of  R  ;  (a)  R  =  1,  (b)  R  =  2,  (c)  R  =5 

o  o  o  o 
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Figure  13.  Total  Cross  Section  of 
Ellipsoid  vs  R  for  0o  =  0°  and 
p  =  2.0 

6.  RESULTS  OF  THE  PROLATE  ELLIPSOID  IN  A  MAGNETIC  FIELD 

We  have  performed  some  scattering  calculation  for  the  case  of  a  prolate  ellip¬ 
soid  with  p  =  2.0  with  the  magnetic  field  aligned  along  the  major  axis  of  the  ellipse 
and  a  gyromagnetic  ratio  of  0.3.  Calculations  of  the  differential  cross  section  for 
both  ordinary  and  extraordinary  rays  with  Rq  =  1,  2  and  5  are  presented  in  Fig¬ 
ures  14  and  15.  The  cross  section  for  the  ordinary  ray  has  a  somewhat  lower 
cross  sec+ion  than  that  for  no  magnetic  field.  The  cross  section  for  the  extraor¬ 
dinary  ray  is  always  higher  than  the  nonmagnetic  field  results  for  RQ  =  1  but  is, 
however,  higher  only  near  forward  scatter  for  other  values  of  RQ. 


7.  COMPARISON  OF  GAUSSIAN  AND  CONDUCTING  ELLIPSOIDS 

To  show  the  effect  of  a  change  from  a  continuous  density  distribution  to  the 
extreme  case  of  a  conducting  body,  we  have  plotted  the  scattering  cross  sections 
at  =  0  and  Rq  =  3  of  the  Gaussian  ellipsoid  (p  =  2)  and  the  perfectly  conduct¬ 
ing  ellipsoid  (Crispin  and  Siegel,  1968/.  The  surface  of  the  conducting  ellipsoid 
was  chosen  to  correspond  to  the  critical  surface  of  the  Gaussian  ellipsoid.  As 
indicated  by  Figure  16  the  cross  section  for  the  Guassian  distribution  is  larger 
for  a  range  of  forward  scattering  angles  but  decreases  so  rapidly  that  it  is 


Figure  14.  Differential  Cross  Section  of  Ordinary  Ray  vs  Scattering  Angle  for 
p  =  2.0,  0ft  =  0°,  and  Several  Values  of  R0;  (a)  R0  =  1,  (b)  RQ  =  2,  (c)  R0  =  5 
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Figure  15.  Differential  Cross  Section  cf  Extraordinary  Ray  vs  Scattering  Angle 
for  p  =  2.0,  en  =  0°,  and  Several  Values  of  R  ;  (a)  R  =  1,  (b)  R  =  2,  (c)  R  =  5 
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Figure  16,  Comparison  of  Conducting  Ellipsoid 
and  Ellipsoidal  Electron  Density  Distribution 


considerably  below  the  cross  section  of  the  conducting  ellipsoid  for  the  range  of 
backward  scattering  angles.  Similar  results  (not  shown  here)  are  obtained  for  the 
case  of  the  elliptic  cylinder  with  the  incident  ray  along  the  major  axis  of  the 
ellipse.  These  results  corroborate  previous  work  for  the  Gaussian  sphere  (Klein 
and  Mabee,  1968),  which  indicates  that  the  replacement  of  a  moderately  overdense 
Gaussian  sphere  by  a  conducting  sphere  of  critical  size  is  a  poor  approximation. 

In  highly  overdense  bodies,  however,  the  density  gradient  in  tl  c  critical  region  is 
quite  high  so  that  here  this  approximation  may  be  a  reasonable  one. 

8.  CONCLUDING  REMARKS 

It  should  now  be  obvious  that  the  shape  of  an  electron  distribution  and  its  ori¬ 
entation  to  an  ooserver  are  very  important  in  determining  the  differential  and  total 
cross  section.  Thus  with  the  figures  available  in  this  report,  it  is  hoped  that  cross 
section  values  of  charge  distributions  can  be  approximated  by  properly  applying 
the  data  presented  and  taking  care  of  end  effects  for  finite  length  cylinders. 


We  have  shown  that  the  magnetic  field  decreases  slightly  the  differential  cross 
section  for  the  ordinary  ray,  but  the  extraordinary  ray  has  a  larger  cross  section 
than  the  same  charge  distribution  when  we  do  not  have  a  magnetic  field.  We  have 
indicated  the  error  involved  by  trying  to  replace  a  charge  distribution  by  a  conduct¬ 
ing  body.  The  characteristic  differential  cross  sections  are  so  different  that  we 
could  only  use  a  conducting  body  for  a  very  small  range  of  scattering  angles.  Al¬ 
though  the  computer  code  in  Appendix  A  was  written  for  a  CDC  6600,  it  should  be 
easily  converted  to  another  system  with  minimal  effort  so  that  individual  cases 
could  be  calculated.  In  addition,  subroutine  DERIV  is  written  for  a  Gaussian  charge 
distribution.  If  another  charge  distribution  were  desired,  only  this  subroutine 
would  be  affected. 


Appendix  A 


Computer  Code  for  Ellipsoidal  Scattering 
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PROGRAM  ELLI®SCINP’JT,0UTOUT,TAI»E5=INPUT,TAPE6=0UTPUT,TAPE39t 

1XAEEZ1 - 

EXTERNAL  0£OIV,0UT> 

COMMON  Xf.3aa>,7.M3Jl»AU,HA.M03).Tr.QUNT _ 

COMMON  0<100>,3CATAM{100> 

nTHFNSTOM  XI  >1-1  tX)f  Yl  IBfmi  .  T  TT  l  £  I M  f  Yl  4BX31  ,X>f?l  tYVI?lTTTTL  fill 
DIMENSION  PRNTm,Y(3»,9ERYf3» 

—DIMENSION  XU0i  f.31.ri  1B1  I3J-»EACXagX4JU _ 

DIMENSION  0IS’.<2),DUNY<2J  ,ANSI(2),RAYC2,?> 


THE  cOLLOHtNG  ARRAYS  DE°EN0  ON  THE  NUMBER  Op 

INIEgEOLATED. MINIS _ 

summon _ ansimixml oa  , amium i .  n roa t  sr  m .  s  tr,i  1 1& a jusnamm am 

COMMON/CONST/  RHO,XO,ZO, THET A, THETAH , YY Y,I I 

MU  (RAY<K>.H.»K:=-U23Z2.aM ORDINARY— R  LI _ / _ 

DATA  (RAY(K,2I,<=1,2)/20H  EXTRAORDINARY  RAY  / 

■D  A-TA.JXTLE/aT.lfl-H _ L _ 

OATA  TITL/34MCR3SS  SECTION  VS.  SCATTERING  ANGLE  / 

DATA  XI  ARI  fYI  A9I  _ L _ 

DATA  XL ADEL/ 3!H  SCATTERING  ANGLE/ 

OA1.A  YL  A-0 EL-Z-Kt-dlM  ?A£-I-.°.AiR  A  MET.ER _ L _ 

“1=3.1415926535896 

^Di-PX/iatl. _ 

TTITL=1 

■ITLsIU - 

XS I TE =11. 0 

■YSIIE-iXl^J - 

I“AGE=  0 


It=  1  pOR  THE  ORDINARY  RAY 

IT  =  2  MR  INF  EXIRAO.RniMARY  RAY 


ITIME  =  THE  NUMOE?  Dp  CASES  TO  DE  PUN. 


DO  501  ITIME=1,4 


READ  IN  HARDNESS,  MAJOR  AND  MINOR  AXIS,  AND  INCLINATION  OR 


_  ELJLIP-SOID _ 


c 

in  cORMaT  (4R10.0> 

_  gsADfg.ini  thttah.yyy _ 

f) 

III  IS  TH F  NHMOPg  DF  DRTPNTA 

tidns  per  c 

ASF. 

f* 

do  61  TTI=1, ? 

11=111 

_LHEIA=T-4£I.A - 

N’=T59 
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IRAG£=I<»AGE*1 

HRITE(6.40)  I°ASE _ 

49  FORMAT  {1H1,12SX,*»AG£*I3) 

0DATE=0ATE(0) _ 

WRITE(6,bB)  OOATE 

&Q  FORMAT  <//,60X.*  THE  PATE  0*  THIS  RUM  IS*AlO./» _ 

WRITER, 20)  RHD»XO, 70, THETA 

— 2JL£0R!liI  tSQX.n  RH3  =  F7. 3,/.58Xil9H  THE  _ 

150X,19H  THE  HIX3R  AXIS  IS,FS.3,7,50X,23H  THE  ROTATIOH  ANGLE  IS,F$ 

_ tsuLrin _ 

WRITE, '6, 110)  {RAY  CK, II) ,K=1,2) 

110  FORHAT  (//7.53X.2A10./7/) _ 

WRITE(6 , 30)  THET AH, YYY 

-J.fl.-^.aRWA,,T._t3JIXi.23-t  _!HI_aAS*lELII_£I£lJL-&MSk E  U*  1x3.0 tl_HlJH-AJ}.TlQlfl 

♦GNETIC  RATIO  OF  F6. 2 ,///) 


0  CALCULATE  ELLI3SOIO 


COST=COSi -THETA* RAO) 

_ SI NT=S IN (-THETA* RAO) _ 

IF{RHO.L£.l.)  RCR=t. 

_ _ IF  (RHQ.LE.  1.  GO  T Q_15 _ 

R0R=30RT ( ALOG(RHO) ) 

_ iHXrJM] _ 

AH7=0. 0 

_  5  00  LI  t-i.lfl) _ 

X1=X0*C0S {FLOAT  {2*TI*RA9)*R0R 

_  Zl=  10.7SlH-(FLJ3A.LL2.yi).*JlAQ)  *RQR 

■(I)=Xl*COST*7l*SINT 

_ lx  q=^x  i*siMl*  Z 1* 1Q.SI _ _ 

AHX-AHAX1{ANX,X{I)) 

_ AH7alHAXl(AM7.?(D) _ 

11  CONTNUE 


“LOT  cLLIPSOIO 


XX(1)=4.0 

_ XX  (1)  =  -4«_0 _ 

YY ( 1) =4, 0 

_ YYi2)_EriuJl _ 

XS  =  1.  0 

_ _ _ 13  =  1*11. _ _ _ 

NO  =  3 

_ IFCIII,FQ.1.AHT.ITIHE.E0.1)  NQ  =  1 _ 

CALL  “LOTV(NO,XX,XS,YY,YS,2,33,-l,ITITL,TITLE,X.ARL,YLA«L, 

_ 1XST7E. VSI7FI _ 

CALL  PL0T^{2,Y,XS,7,YS,181,33,1,ITITL,TITLE,XLAREL,YLA«EL, 

_ 1XSIZE,  YSI7E) _ _ _ _ _ 

C 

_S _ .CAI.CLJIATE  THPA-iT  PAOAHFTFRS. _ 


IF(TH£TA.EO.O.fl.OR.  THETA. EO. 90. 01  GO  TO  25 

_ T?inH=45-TFTXiHEJ47?.> _ 

AXG=130.-TKETA 

_ AHX=(7(THDH)-7(TH0H-1H  7(X(TNOH)-X(TNOH-1)  I _ 

00  31  I=l»50 

_ mQ5ilMS±2*  TEL3A.I  i  DJ  12m _ _ _ 

Axi=<ra*iDj-z<iNO-i»>/(xaNn>-xciNO-i»i 

_ trUMI  .CT.AHX)  TMOHsTHO _ 

IF(AX1.GT.A-(X)  AHX=A4i 

_ 3i  nattmiuE _ 

DATA  FA3T0R70. 010,0. 025,0.050,0.10,0.25, 0.50, 1.0, 2. 0,5.0/ 
2B _ 30.41 -I  =  U1M _ 

O(T) =0.0 

_ 11  CONTINUE _ 

ZER0=7 (INOM)  +0.  3*A*(7*siX(-2.  0*THETA*RAD1 

_ 1F-11M£.T-A^£Q»1UXI..Q.2»THET  A.  F3..90. 0)  7FPn=0..Q _ 

AH7Z=A*T 

_ T.F.:AH7«l£.in*SlS.145. »«AQ) .A NO. XO/Zfl. RE. 3).  AH7sAM7»0. 3 _ 

00  71  1=1,49 

_ aiaaaZ£go _ 

J=8 

_ IE.a.L£.2S)  J=7 _ 

IE (I. LE. 22)  J=3 

_ IFd.LE^lA)  -J?5 _ 

IFd.LE.15)  J=<* 

- IF.ir.LE.17.l-J.a3 - 

IF(T.LE.fl)  J=? 

_ IF  (I. I.  E..4).-Jsl _ _ 

oif=fa';tor(J)*a>(7 

_ IB. 1 0.T.F  .jGT  at.  2-51— DIF  =fl.  23 _ - 

IF (0(51-1 l.LE.'HZZ.ANO. OIF. 5T. 0.15)  OIF=0.15 

_ IE  (01 51-r.I AHO.ai Sir IU I.  AMZ.7U1.  15 _ 

1. AND, THETA. EO. 30 • )  TIC=0.525 

_ 9 15jUJJ.-T.fl  (.51  -X)  t  3.1  F _ 

0(50+1) =9(49+1) -OIF 

_ LL.aQN.tXHU£ _ 

J=  0 

_ on  at  Lai»U3 - 

IF (0(1 ) , EO.O. 0)  GO  TO  81 

j=.m 

_ iijuaui _ 

81  CONTINUE 


I xr»=  MUMOER  Or  IMO10T  oftojMFTERS. 


I0O=J 

_ M3UE(6..SQl-13(l-l.l5l,IflPl _ 

50  fo°8AT  (77.27H  THE  IMOfl:!  0AR4»ETE°S  A*E, / , 10 ( 10F10. 3,/) ) 


C 


START  RAY  TFAdNG 


|c*> 


PRHT(3)=  INITIAL  STEP  SIZE.  IF  IT  IS  NEGATIVE,  EACH  STEP 
IS  NOT  PRINT 


PRNTC4) =  NAXINUN  AX3ULAR  CHANGE  ®ER  STEP. 


pRNT;3)=-3.25 


I 

_ I 

- (PRNT ( 
°A0£= 1° 

3I.LE.C.0)  3 
AGE+1 

0  TO  35 

WRITEI6, 

40)  IPA5E 

1XSI7E.YSI7EI 

3  A T AN ( II = AL°HA ( I30UNT) /PAO 


IF(SCATAN;II.LT.0.3)  SCATANIII  =360. *-SCA TAN  III 


0  TO  SI 


53  IF (S3ATAN I I) < 

r. 


51  IP(SCATAN(I*U.GT.SCATANII))  GO  TO  92 
»1=TN“1-1 


00  101  J=I  , IH31 


ffTWraTTciTIa 


pi  j)=p{j*D 


IF(I.GT.INPl-l)  no  TO  5 
co  to  «;* 


32  IF (IE, NE. 0)  GO  TO  91 


91  30NTINME 


OATC  ANGI/180. ,180.1/ 


WRITE  16 , 10  0)  OISPL.'l) 


iiJMirif  KiMwivii.iini 


iRniUidtvi  *kaww**>  i  wirx+\ i;»i 
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WRITE  TN®‘JT  PARAMETERS,  IMPACT  ®ARAMETCPS,  AND  SCATTERING 

ANOl  r<?  (IN  TA®-?- 

WRITE:?, 70 1  RH0T*nt7ntTMFTAf TWFTAHfYYY, nnATF.TMP 

70 

FORMAT (jE10. 5, AlO, 110) 

WRTTF(7.A(J)  (BUI  ,  SO  AT  AN  (D  ,  r=l,TN®» 

r* 

80 

FORMAT  (8F13.1) 

•% 

r 

CONVERT  SCATTERING  ANGLE  To  RADIANS  AND  INTERPOLATE  B  VS. 
SCATTERING  ANG-E. 

On  131  I =1 

*>(I)=B(I)-OIS».(l) 

SrATAN(TI=SCATAN(T)»PAO _ 

131 

CONTINUE 

OT=?.fl*oT/3ftQ. 

ANOINT (1)=0. 001/57. +OT 

00  ??1  T~? ,N3 _ 

??!_ 

ANOINT <1 >= ANOINT (I-l»*OT 

OONTTNIIC  _  _  _  -  _  - 

CALL  INT°( SO ATAN,B, IN3, ANOINT,  BINT, ODD A  ,N°,JL,J-i> 

00  Ifcl  T s.ll  , .IN 

ANGINT(I-JL«-1I=  ANOINT  (II /RAO 

bt*jt;t-ii  ♦i»=otnt:t»  _  ...  . 

-1-41- 

OOOA<I-JLlil=030A(I» 

.  CONI.THUE - 

C  “LOT  IMPACT  PARAMETERS  VS,  SCATTERING  ANGLE. 


XX(1)=0.0 

-XX.I2I.g36a*Jl _ 

VV(H=-3.0 

yy (?) =3. 0 _ _ _ 

SX  =  (fO. 

_ ST-=JL».fl _ 

CALL  “LOT Vt 3, XX, 40. ,VY,1. 0 , 2 , 33, - 1, ITITL , T ITLE ,XL A«EL , YL ABE L, 

1XSIZE.YSLZ.E) _ 

CALL  “LOTV(2,AXGINT,40.,BtNT,1.0,N“,33, 1,1 TITL, TITLE, XL ABE. ,YLA  BEL 


_C _ SHOO THE.  OERIVA  U.UES _ 

i) 

_ rFUTRaa - 

N“  1  =  N°- 1 

-  no-zoi-MS'i-ivi - 

Ic (NSM, EQ. 1)  GO  TO  45 

_ 00_taL- _ _  _ _ 

IJ=IFILTP 

- lF.(j.LE.lEii,n; — i — I  j=-  - 

IF(JHFILTR.GE.NP)  IJ=N“-J 
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SO»»OA<J>=g.O 

-03-131  JJrJLiIJ 


soboai j>  =soboa ;  j)  obda  ;  j*jjj  oboa  i  j-jj> 

191  CONTINUE _ 


S090AC J)=SOBDAtJ)*OBOA<J) 

_ SQ90A(Jt=S09QA(  UZ£LvMi2^LLty_ 

m  continue 

_ 03.17.: .  _ 


090A<J)=S090A(J» 
in  CONTINUE 


IFINSN.EO. 2>  S3  TO  201 


CALCULATE  CROSS  SECTION 


45  SUM=3.3 

_ D3  .121  Irlt’lg- 


SIGL  tl )  =4»  0*<*I*BINTlI)  *090 A ( I)  /SI  NC  ANGINTI I)  *RA3) 
121  CONTINUE _ _ _ 


00  111  I=1,NP 

■SO =1.0 _ 


IFCt.EQ.l.OR.I.EO.NO)  SC=0.5 
SU*=A9S<SIGL(I)  ) *SC»SJ* . . 


SISLCIMALOSlO.'APSCSIGLam 

.111  GMU.NUE _ 


SUN=SUM+0.5*(*ANGINT<1!  »»<0«0»(1)  * (- ANGINT  1 1) )  ♦10.**SIGH1>  » 
SU.N=SUN*0.5*(1SQ.-ANGIHT(NP)>»(OO0AtRP>*(ANGINTtNP)-a&3.>n0»** 
1SIGL JN°) I 

_SU?1.-mHJL0AD _ 


WRITE<6,120>  SJi 

.121 . fOSJiAT  ( _/// .  30< « *  TOTAL  CPOSS  SECTION  =»E9.3«8H  (KM*»2>» 

IPAGE=I“AGE*1 

_ w*.n.Et6.<tQ)  ijv,e _ 

L=NP/100*1 

_ 0Q_  L^SEljtJ. _ 

LL={R-1>*100U 


IF(K.EO.L)  LU=N° 

-H.RIT.E16.i90l  ...  .  - 


90  FORMAT  J 2!*  SCATTERING  ANGLE 
_ LG_SI£NA1j3JUJ _ 


IMPACT  PARAMETER  OERIVATIrfE 


-( 


WRITE <6, 170)  ( ANSINT (IT , RINTtll , OBOA ( I ) ,SIGL 1 1> , I  =  LL ,LU» 

_ 1ZJL-EQLPMI1 121aiU-E  /.«Z»13  X^ES.  UliiT  LmJ*  jSXif  6*Jjl6XU _ 

IPAGE=I°AGE*1 

_ MlIIElfultm  IaftiE _ 


151 


CONTINUE 

jixiii-=a^iL 


xx;2)=Too.o 

..miL-JuJ _ 


Yv ( 2) =-4»  0 
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SUBROUTINE  D£RI/<YY,OERY,NEG» 


C  THIS  SUBROUTINE  CALCULATES  THE  DERIVATIVES  OF  AN  ELLIPSOI3AL 

11 _ CHARGE  DISTRIBUTION  FOR  THE  RUNGE-KUTTA  HETHOO. _ 

Hfc 

CQNHON/CONST/  RHQ.XQ tZO« THETA.  THETAH.Y. I _ 

OINENSION  H(  2)  .  3  MO  A (2I»YV<31 »OERY  <31 ,0UDX<2» tOUQ?<2) 

_ "1=3. 1415Q265 _ 

RADrPI/180. 

_ NE6=3 _ 

COST=COS<THETA*RAD) 

_ SINT=SIM<THETA»RAO) _ 

CHI=«?H0*EXP(-(<YY<1)  *CnST+YY(21*SINTl/X0)**2 

_ 1-t t-Yv;i)»SINT»YY;2)»C0STl/7Ql«»21 _ 

YL=Y*OOS< (-THETAH’RAOl  +YV<3)  1 

YTgV»SIH((-THETAH»RAQ)»YY(3)l _ 

HI  =  <YT»*41/4.MYL**21M<1.-CHI1**21 

—  if;wi.lt.o.o)  ia  to  5 _ 

HI  =  SDRTIH11 

—  Hti;=l.-CHI-(YT»»?!/2.+Wl _ 

H< 21=1. -CHI- <YT**21 /2.-H1 

_ ahu=i.-chi»;i.-chi)/h;d _ 

IFUHU.LT. 0.0!  SO  TO  S 

iFm.ea.o.)  quo  a  =  o.o _ 

IF (HI. ED. 0.1  S3  TO  35 

OWOA:i)=VL»vT»;-i.  +  ;vT»»2-?.»a.-CHi)»»2)/{2.»Hiu _ 

OtnA<2)=YL*YT*(-l.-<YT*»2-2.*< l.-CMI! ** 21 / ( 2. *Wll 1 

_ QUO A=CHI» < l, -CHI)  »OWOA  <1!  /  <9. »S3RTUHU1  »  <W  (!)»»?! ) _ 

35  CONTINUE 

_ ochiox=-2.»chi» : ;yy : n  *oost»yy ;2i  »sinti »cost/ ;xo*»2) _ 

l-(-SINT*YY(lH-YY<2)  *COST!*SINT/(TO**2)I 

_ acHiRTsA  »chi;  u  m  n  ^aasmmiJLSiim  »^int/ ;  win _ 

1*  <  — YY ( 1) *SINT*yy<2) *C0STJ *C0ST/<70**2>> 

_ IF.;h1«.£1.  0.Q1..  »Q  TO  «t5 _ _____ _ _ 

OJOX'l)=nCHIOX*(2.*-;HI-t.-CHI*  (l.-CHH*(l,  +  (l.-CHI)*(YL**2)  / 

ouox<2)=o:hiox*<2.*chi-i.-2mi*;i.-chi)*{i.-;i.-:hi)*;yl**2) / 

_ lHl)/HIHl/<2.*a3RT(AHUl  *H(H1 _ 

nun7<ii=n,CHlo7*<2.*CHi-i.-DHl*(i.-0Hi)*<i.^<i.-3Hi)*<YL**2)  r 

LHl)/H<T))/:?.*50RT;aHU)»w;i)! _ _ _ 

DUO7 (2)  =  DCHI07*  <?.*0HI-1.-0HIM1.-CHI1*<1.-<1.-:HI)  *<YL**2>  / 

_ 1N1)/H<H1/<2.»30RT<AH'J1»H(I)1 _ _ _ _ 

GO  TO  55 

_ 45  nunx(D=-o.CHnx/  ;7.»sort  uhuii _ 

OUOX<?)  =-0CHI")X/<2,  *SORT  <  AHU) ) 

OUO7<2>=-0CHn7M2.  *SOdt<AHU>> 

_ 55  YY3  =  YY:3» _ _ _ _ _ _ 

OXOT=(SDRT(AHLM  *''0S(YY3)  +OUOA* SIN < YY 3!) /AMU 

_ 07nT=<50RTf  AM.J)  »STN<YY3)-0U0A»0nS(YV3)l  / AMU _ 

OAOT=(0OS(YY3)*OUO7<I>-SIN<YY3)*OUQX<Il )/ AMU 


SO  TO  15 

_ 5_bl£IUrJ _ 

SO  TO  25 

_ is.  ai£iML=  nxni 

0£PY<?) =Q7DT 

_ O  F  O  Y (3 >  =OAQT 

25  CONTINUE 


return 

END 
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subroutine  out:»{xx,y,o£bv,ihlf,  prnt,neg> 


'  this  SUBROUTINE  PRINTS  THE  9 "SUL TS  OF  the  RUNGE-KUTTA  INTEGRAT  IDN 

1- _ ANO  OTTEgilMES  WHEN  WE  HAVE  FOLLOWED  A  RAT  FAR  ENOUGH. _ 

_ DIMENSION  v:3).0ERY(3) .°RMT(5) _ 

T- <0£RY{3) .LT.l. OE-5.ANO.XX.GT.O. 0)  fpmt (5) =*»1#0 

_ 15 i «*.?! <J»tl.T...QJfl>-SQ-T9.  5 _ 

WRITE(6,10)  XX,  (V(Jt ,J=1,3) , (OERY (J) , J= 1, 3) , IHLF, MEG 

10  FORMAT  (7E17.8,2TS> _ 

5  RETURN 

_  _IV  _ _ _ _ 


SU^OUTINE  R<SS;PRMT,VV.O"RY,Mf)IM,THLF,OERIV,OUT'»> 


THIS  SUBROUTINE  ’ERFORWS  THE  RUNGE-KUTTA  INTEGRATION, 


DIMENSION  PRMT.-l)  ,TY  ;i)  , OERV ;i) , X 1 {3 03)  ,t;330>  *  ALPHA  1300)  ,<  H,V» 

_ 0 IMENSJLQN.  JCX  130  OJUQA  01141 _ 

REAL  < 

_ OO.MHON  XI ,7 . AURA  A ,J _ 

OOX^ON/OONST/  RHO,XO, 70, THETA, THETAH,YYY, III 


SET  ijo  THE  INITIAL  CONDITIONS. 


°T=3. 14150265 

RAD^I/iaQ.. _ _  _ 

THET  =  THEMH 

IC(AQS(T HI TAH)«3  T»9Q»  JlIHET =.ia.a^=JiaSXDlEIAlU _ 

T^£T  =  THET*Pn 

rio=o  . . . . . . __ . . . 

,J  =  0 

oion  =  a.o  .  ... _ 

IHLC=3 

HsAOSC^RNTO)) 

_ IRS=.o.0.-M.T  (41.  _ _ _ 

PRNT (5) =0. 3 

_ 3a  J=J+L _ 

Ic (J.GE. TOO)  H3IT£(5,10) 

_ J.0  FORMAT  l//,2aX.*IH:  NUMBER  QP  STE°S  HAS  FX  SFFDFn  THE  DT.MENE  T.OM*  I 

IP«J,5E*300)  GO  TO  55 

xi(j)srrd)  _  _ _ _ _ _  _ — _ , _ 

7< J) =YV(?) 

_ -Alp.HALJ1  =  TTOL _ _ _ 

xx  ;j)sv 

_ 45  11*0  _ .  .  _ , _ - _ 


GALBUHTE.  THE  DERIVATIVES,  AND  .  14  E EX  .THAT  JtQ_£RRQ3_J lAS-.OSlUlRUI . 

IN  0E°':y. 


5  GAUL  OERTV(YY,0£RY,NEG) 

GALL  QUT°(X,TY^0ERY,HLc,:»oMT,NEG) _ 

IC(P°MT (5) ,NE«  0. 0)  GO  TO  55 

.  _  .  -  .  n=ii+i _  _ _ _ 

IL ( I  AOS ( NEG) • GT, 0)  GO  to  15 


P  <r> 


STORE  CONSTANTS  ANO  CALCULATE  THE  NEXT  SE  OF  VALUES  AT  MHI^ 

■er 


K( I .  II) =OERY(I) *H 


OAOT(II) =DERY( 3) 


1^(11. FI. 3) 
yr<i)=xi;j)*K 


YY(2)=?(J)+K12»II)/C 
YY(3)  =  ALPHA(J>*K<3.II)  /C 


X=XX< J) +H/C 
5 


THE  DERIVATIVES. 


15  H=H/2. 


Ip  ( IHL c. GT .40)  50  TO  55 

-I  r.  1=  I- 


x=xx;j> 


VV(2)=X(J) 


GO  TO  45 


CALCULATE  THE  VALUE  OF  TH"  FUNCTION  USING  THE  FOUR  PREVIOUSLY 
CALCULATED  p 


0ELZ=;k,-2,1»*E.*  ;<;2,2>*KJ2,3)  )♦<  12,4) >  V6. 

0£LA=(K(3.1) +2. * (K (3.2)+K ( 3 . 3) ) +K ( 3.4))/6. _ 


DAOT2=AOS<<OA0T<1>«-OAOT<2M-OAOT(3)*OAPT«4)  >  V4.) 

~ CHECK  that  THE  ANGULAR  CHANGE  IS  NOT  GREA ter  THAN  OESITEO, 


I c (AOS (DEL A) .ST.  ERR)  GO  TO  15 


STORE.  ACCEPTAPLE  VALUES  OF  THE  INTEGRATION. 


YY(l)=Xl(J)OE.X 

XJ 12)  =  7A±)*3&  7 _ 

yv  ;3)  =  al»ha;j)  *dela 

-JL=xx.ui+a _ 


IF(RHO.LT.l.fJ)  SO  T0  11 

_ IElIMQ..£ja.-U_-S3  -IO-U _ _ _ _ _ 

RCR=SORT( ALOG(RHO)) 

_ LFA1XJS1.*±Z±X1111*1Z.  -RC.R  _ 


IF(APS(YY(3)).-T,APS(THET))  ER  =(A°S(YV(3))”APS(THET))*2. 
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r 

°HASE  ANGLE  HURT  RE  WITHIN  0.?5  OEGOEE*  OF  THE 

MAGNETIC  FIELO. 

IF(4RS(ER).GT.R4r'*.20)  GO  To  u 

YYf3>=  YYnHSf.Uffa  .YY(31  1 

INO=l 

0 

CHECK  IF  THE  INTERVAL  OF  In^-GRATION  CAN  RE  INOREASEO. 

11  Ic(ARS(OELA).LT.O.ai.ANO.OAOT?.Lf.OAOTl)  GO  TO 

E5 

OAOT1 =DA0T2 

GO  to  3*- 

«R  OAOT1 rOAOT? 

IF(H.GT,  AOSI’R'IT  (3) >  >  GO  TO  33 
_ a=a.*2i. 

IHLF=IHLF-1 

_ 110.  TO  33. _ c _ 

5 3  RETURN 

_ Em _ . 
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